Single and two-particle motion of heavy particles in turbulence 
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We study motion of small particles in turbulence when the particle relaxation time falls in the 
range of inertial time-scales of the flow. Due to inertia, particles drift relative to the fluid. We show 
that the drift velocity is close to the Lagrangian velocity increments of turbulence at the particle 
relaxation time. We demonstrate that the collective drift of two close particles makes them see 
local velocity increments fluctuate fast and we introduce the corresponding Langevin description for 
separation dynamics. This allows to describe the behavior of the Lyapunov exponent and give the 
analogue of Richardson's law for separation above viscous scale. 
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Motion of small particles in a fluid, due to random 
molecular forces, is the subject of the classical theory of 
Brownian motion. Scale separation between the particle 
relaxation time and time-scales of the forces allows to in- 
troduce an effective Langevin description of the driving 
force as white noise in time In contrast, here we con- 
sider the situation where the random driving force origi- 
nates not from the microscopic motions, but rather from 
the macroscopic turbulent motion of the surrounding 
fluid 0] • The limit of the particle relaxation time much 
larger than the characteristic time-scales of turbulence 
(very heavy particles) can be described as in the Brown- 
ian motion case Q. In the opposite limit, when particle 
relaxation time is much smaller than the characteristic 
time-scales of turbulence (very large friction), particles 
follow the flow closely and the two-particle dispersion - 
of interest to us here, - is approximately the same as 
for fluid particles. In this Letter we study the interme- 
diate case of heavy particles, where the relaxation time 
falls in the range of flow time-scales corresponding to the 
inertial interval of turbulence. This precludes Langevin 
description for the single-particle motion. However, for 
two particles, because of their collective drift relative to 
the fluid, velocity increments determining the separation 
do vary fast. This allows to introduce effective Langevin 
description for the dynamics of separations. The descrip- 
tion enables us to find several results on particle behavior 
in turbulence valid beyond Kolmogorov theory. 

Behavior of small inertial particles in turbulence has 
received much attention lately 0, H, @, 0, H S 53, E 



3 B Q E3, EE [HI- This Problem has marry appli- 
cations including rain formation in clouds 0, ocean 
physics [6} and engineering [7j. Theoretical study of the 
problem mostly involved modeling turbulence by a white 
noise in time, Gaussian velocity field, known as Kraich- 
nan model. Even in that case theoretical study is rather 
difficult, analytic results were mainly obtained for the 
one-dimensional case [1, 0]. The limit of heavy particles, 
considered here for turbulence, was studied numerically 
for Kraichnan model in 0, E3] ■ For turbulence, numer- 
ical studies of intermediate regime of moderately heavy 
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particles were performed in 

We consider the motion of a small spherical particle 
in an incompressible, statistically steady, turbulent flow 
u(x, t). We assume that the drag force acting on the par- 
ticle obeys Stokes' law. Designating the particle position 
and velocity by x(t) and v(t), Newton's law reads 



v = — 



v — u(x(t), t) 



(1) 



Here r — (2/9)(po/ p)(a 2 /v), where po and o are the par- 
ticle density and radius, while p and v are the fluid's 
density and kinematic viscosity. We briefly review rele- 
vant properties of u(x,t), see e.g. [2| for details. Velocity 
field, excited at the integral scale L, fluctuates in a wide 
(inertial) range of spatial scales rj -C I -C L. The char- 
acteristic velocity it/ of fluctuations at a scale I is related 
to the temporal scale ti by uiU/l ~ 1. At the viscous 
scale rj we have t v ~ rj 2 /v. For moderate Reynolds num- 
bers one can use Kolmogorov theory (below K41) that 
gives m ~ (el) 1 / 3 , U ~ e -V3/2/3 and ^ _ (j/S/ejV*, 
where e is the mean energy injection rate. Eq. ([1]) is 
valid if r] 3> a and inertia-induced particle drift relative 
to the flow, described by w(t) = v(t) — u(x(t),t), has 
small Reynolds number wa/v [l6j]. We shall consider 
t n <C t <C which implies that particles are heavy, 
paj p ~ (rj 2 /cl 2 )(t /t.^) ^> 1, justifying the neglect of such 
effects as added mass in Eq. ([1]) [16j. Beyond K41, quan- 
tities like 77 and w have strong spatiotemporal fluctua- 
tions, and in that setup we will refer to their local (in 
space and time) values on statistically relevant events. 

We first consider the drift velocity w. From Eq. (JTJ) , 
in the steady state, w = f_ uexp(t/r)dt/r, where 
5fu = u(x(t), t) — u(0, 0) is the turbulent velocity differ- 
ence in particle frame and we chose t = and x(0) = 0. 
The difference 6f u is analogous to the Lagrangian differ- 
ence 5^u = u(q(t),t) — u(0, 0) where q(t) is a fluid parti- 
cle trajectory obeying q = u(q(t),t) and q(0) = 0. Just 
as S^u, the increment Sf u should be, on a rough scale, a 
non-decreasing function of \t\ , growing at most as a power 
law. Then the integral for w yields the order of magni- 



tude estimate 



and we also see that the charac- 
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teristic time of variations of w is r. Let us show that from 
w ~ S^u one can pass to w ~ <5^u, that is <5:f it ~ 
We write the telescopic sum = S^u + Sv, where 

6v = u(x(t),i) — u(q(i), t) is the velocity difference of the 
separating fluid and inertial particles. The separation is 
due to the combined effect of the inertial drift and the 
explosive separation of trajectories in the inertial range 
[l9( |, so that Sv < max[w, 5fU], where S^u is the relative 
velocity of explosively separating particles 0. Thus, at 
t such that w < §f u, one has 6v < 6^u and, then, from 
the telescopic sum, 8fu^ 6fu. This, combined with 
w ~ <5^w, allows to show iu ~ S^u. Since Jfu roughly 
grows with t at < t < tx and = <5^ W < 10 ~ ^t L u 
(the latter is just an upper bound for w), there exists i» 
such that w ~ 6f u. At this £*, from the previous con- 
siderations one has <5^u ~ 6p,u ~ iu ~ <5^m. Finally u 
roughly increases monotonically for timcscalcs less than 
L/w > L/8t L u ~ tt, so we conclude that t* ~ r, whence 
it) ~ <5^w. This means that particles follow only the flow 
fluctuations with time-scales larger than r. Though ex- 
pectable, the result uses specific properties of turbulence 
and does not hold for any u(x,t), where the relation 
between the spatial and temporal fluctuations is differ- 
ent. For the particle velocity, from w <C u, we have 
v w u(x(t),t). Note that w is accessible experimentally 
through acceleration a = v = —w/t. Also note that 
w ~ 8^u applies at any r: at r > t^ it is trivial, while at 
t trq one has w w — r[<9 f it + (it • V)tt] ~ 6^u (IBj . 

The local equality w ~ 8^u suggests that the time- 
averages along the particle trajectory x(t) satisfy 

r>»> = („,») ~ ((#„)«) ~ ul(r/t L )^. (2) 

The anomalous exponents 7„ could differ from their 
counterparts for the more usual ((8^u) n ) q , where the 
averaging is along the fluid particle trajectory 0, lioj ]. 
- inertial particles tend to concentrate preferentially 
in specific regions of the flow [4j. Yet the difference 
seems unimportant: introducing l T by ti T ~ r, one has 
{{5r u T) ~ {[Su(Lr)] n ), where 5«(J2) = u(x + R)-u(x), 
and ([<5ii(Z T )] n ) allows additional spatial averaging over 
l T vicinity of x(t), beyond which the preferential concen- 
tration is expected to be small [13|. In K41, similarity of 
properties of w and 5^u follows by dimensional analysis: 
one has w ~ ^/er and (w n ) ~ (er)™/ 2 like for (5^u (20 |. 

We now consider two-particle motion. We assume the 
particle separation R — x'~x much larger than radius a, 
so that hydrodynamic interactions between particles are 
negligible. Then each particle satisfies Eq. ([T]), producing 

R + R/t = 5u(R)/t. (3) 

At R 3> l T particle dispersion, driven by 8u(R) ~ ur, is 
determined by turbulent fluctuations slower than r, and 
the separation is like for fluid particles: R — v' — v = 
8u(R)+w' — w w Su(R), by w ~ u; T <C Ur. In contrast, 
at R <C i r dispersion laws peculiar for inertial particles 



hold (in K41 Z T - e^V 3 / 2 - //(r/t,,) 3 / 2 ). We first treat 
i?(i) <C ?7 where Sui(R) w RjV jUi(x(t) ,t) and 

R + R/t = (R-V)u/t. (4) 

The main characteristics of Eq. ((4]) is the Lyapunov expo- 
nent Ai = limi_oo ln[p(i)/p(0)]/f = (p/f>) describing the 
exponential growth of the distance p = (R, tR) between 
two infinitesimally close trajectories in the phase space 
of Eq. (fTJ). Here Ai also describes the growth of R and R 
separately, so it can be observed via spatial trajectories 
of close particles with small relative velocity. 

We assume that the time-average (p/p) can be found 
by averaging over the statistics of turbulence: Ai — 
limt_ >00 (p(i)/|)(t)) u . It is useful to consider first the 
Kraichnan model where V ' jUi in Eq. ((4]) is modeled 
by white noise &ij obeying (uij (t)a mn (t')) = D5(t — 
t') [(d + l)5 lm Sj n - 5ijS mn - 8 in 5 m j] where d is the space 
dimension [6;]. Passing to dimcnsionless time s = Dt, one 
finds (Dt)R + R = a'{s)R, where (a' t3 (si)a' mn (s 2 )) = 
S{si - s 2 ) [{d+ l)5i m 5jn - Sij5 mn - 5inS m j]. At Dt <c 1 
one may drop (Dt)R. The resulting equation is the 
same as for separation of fluid particles so Ai w A{', 
where A{ ~ D is the Lyapunov exponent of fluid parti- 
cles [lj|. On the other hand, using dimensionless time 
a' = D^H/t 2 / 3 , one finds R + R/iDr) 1 / 3 = a'(s')R, cf. 
[n|. At (Dt) 1 / 3 > 1 one may drop ^/(Dt) 1 / 3 . Then, 
by dimensional analysis, Ai ~ D 1 / 3, /t 2 / 3, . We observe 
that at Air <C 1 one can drop the first, inertial, term 
in R + R/t = <jR/t, while at Air ^> 1, the friction 
term can be dropped: the characteristic time of varia- 
tions of R is Aj" 1 so the ratio of R to R/t is estimated 
as Air. The region of small inertia Air ~ Dt <C 1 is sepa- 
rated from the region of large inertia Air ~ (Dt) 1 ' 3 ^> 1 
by the long crossover region Dt > (Dt) 1 / 3 ~ 1 where 
R ~ R/t and Air ~ 1. This explains the numeri- 
cal results of 0, [l(|. For example, in d = 2 one may 
write Air = [Dt^X^Dt)- 1 / 3 ], where Ai(e) = QRz) 
for complex dynamics i = — z 2 — ez + 7 [6(. Here un- 
correlated noises ji obey ("fi(t)ji(t')) = 8(t — t') and 
(l2(t)l2(t')) = 36 (t — t'). The crossover region is de- 
scribed by slow, order unity, variation of Ai(e) from 
A(l) w 0.5 to A = Ai(0) pb 2 3], and Ai w X°D 1 / 3 /t 2 / 3 
at (Dt) 1 / 3 > 1. The time-scale beyond which (p(t)/p(t)) 
relaxes to its steady-state value Ai, forgetting the initial 
conditions (for d = 2, the relaxation time of z), can be 
estimated as A^ 1 . Indeed, Aj -1 is the only time-scale both 
at Dt < 1 and (Dt) 1 / 3 > 1, while at (Dt) 1 / 3 - 1 all co- 
efficients in R + R/ (Dt) 1 / 3 = <r'(s')R are of order unity 
so again t 2 / 3 / D 1 / 3 ~ A^ 1 is the only possible time-scale. 

Let us now consider the turbulent velocity gradient 
seen by the particle, V jUi(x(t),t) in Eq. It is deter- 
mined by fluctuations at the viscous scale, Vu ~ u v /rj, 
while its time variation can be inferred from 

^-Vu(x(t), t) = [(d t + u ■ V) + w ■ V] Vw. (5) 
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At t 3> t. n , the drift derivative (w ■ V)Vw ~ Vu(w)/ry) 
dominates the substantial derivative (d t + u ■ V) Vu ~ 
\7u(u v /r]). Thus \7 jUi(x(t), t) varies at time-scale r\jw - 
during this time particle deviates from the carrying flow 
by the spatial scale of variations of velocity gradient, 77. 

Below we study Ai as a function of the Stokes num- 
ber St = A* r. Here A'" 1 " 6 is the Lyapunov exponent of 
fluid particles in turbulence. K41 dimensional analysis 
gives Aj nrh ~ \Jt/v, while it does not fix Ai due to the 
additional time-scale r. At physically relevant Reynolds 
numbers X[ urb ~ y/ej^ is valid [H, indicating that X[ urb 
is determined by weakly intermittent events, and imply- 
ing the same for Ai at least for small St where Ai « Af urb . 
We shall assume that Ai is determined by weakly inter- 
mittent events for any St and use K41 for order of mag- 
nitude estimates, justifying this later (e. g. t v and 77 will 
refer to their K41 values). The domain of applicability of 
Ai is St 1 - at these St the drift contribution in 

Eq. (O is negligible, while Air <C 1 allows to neglect the 
inertial term in Eq. ([!]). We now consider St 3> 1. The 
idea of the analysis is to make the natural assumption 
Ai <C X[ urb , easily verifiable in the end of the calculation, 
and to use that the timescale A^ 1 of variations of R is 
much larger than the time-scale of variations of Vu in 
Eq. 0}: A7 1 > [\{ urb ]- 1 > i]/w. Then the effect of Vtt 
on R can be represented by a white noise. 

The first step is the derivation of the Lyapunov expo- 
nent Ai (a) for the auxiliary constant-drift problem i?j + 
Ri/r = RjVjUi(x a (t),t)/T, where x a = u{x a (t),t) + a. 
Having in mind the application, the constant vector a 
is assumed to have the same characteristic value -y/er 
as w, so that the correlation time r c of 'Vu(x a (t),t) is 
77/a ~ l/A' ur6 VSt. Averaging Eq. (gj) over time At sat- 
isfying t c <C At <C AT/, we find i?, + Ri/r = RjX7 jUi/r, 
where V j Mi = J 4 ' +A * V jUi(x a (t'),t / ) dt' / At. We observe 
that VjUi is a Gaussian process with zero mean and pair 
correlation which - due to stationarity, spatial homogene- 
ity of small-scale turbulence and incompressibility, - is 
determined by Fij mn (a) = Jdt (VjU^O, 0)V„u m (q(t) + 
at,t)}. Then \7jUi(t) is statistically equivalent to 

Jt +At crij(o,,t')dt' I At, where <Tij(a,t) is a white noise: 

(aij(a,t)) = 0, (aij(a,t)a mn (a,t')) = S(t' '-t)F ijmn (a) . 

Dropping the auxiliary time-averaging we conclude that 
Ai (a) can be found from the anisotropic Kraichnan model 
R + R/t = oR/t. For 77/a <C t v the time variation of 
\7 jUi(x a (t), t) is determined by the drift and one can 
simplify Fij mn (a) J dt(V jUi(0)V n u m (at)) . The aris- 
ing degeneracy a n F ijmn = J dtd t (V jUi(0)u m {at)) = 
allows to set 0^3 = 0, where we chose z— axis parallel 
to a. As a result, r = [R\,R^) satisfies closed equa- 
tion r + r/r = err/r where a is a 2 x 2 matrix with 
<Sij = aij, while i? 3 = <7 3i ri. Using (VjUj(0)V„u m (r)) = 
V j V„S , im (r)/2, where Sy(r) = (Mr) - Ui(0)][uj(r) - 
Uj(0)]), one can express Fij mn with the help of sec- 



ond order structure function of turbulence S^M = 
(([u(r) — u(0)] ■ r/r) 2 ). After straightforward calcula- 
tion one finds that a obeys the statistics of 2d Kraichnan 
model with D = D{a) = a" 1 / Q °° S 2 (r)/ (2r 2 ) dr. Thus 
Ai(a) = D(a)V 3 A\[(£(a)r)-V 3 ]/T 2 / 3 , for 77/a < t v or 
\/St > 1. In general, Xx{a)r ~ (D(a)r) 1 / 3 at St > 1. 

We now show Ai ~ (Xi(w)}, where the average is over 
the single-time statistics of w. If we let a above be- 
come slow function of time, then p{t) grows with local 
Lyapunov exponent, (p(t)/p(t)} w Ai[a(t)], provided the 
characteristic time A 1 " 1 (a) of the relaxation of (p(t)/p(t)) 
to its local value (see the discussion on Kraichnan model) 
is much smaller than the characteristic time of variations 
of a(t). Next, the drift velocity w(t) is determined by 
velocity fluctuations with ti > r and thus its fixation 
influences weakly the statistics of velocity gradients de- 
termined by fluctuations with time-scale t v <C r. Thus, 
(p(t)/p(t)) for Eq. ([4J can be found first averaging over 
the gradients at fixed w(t) and then averaging over w(t). 
The first averaging gives Ai [?«(£)] provided Ar/ 1 [w(t)] is 
much smaller than the characteristic time r of variations 
of w(t). However, under the latter condition Kraichnan 
model gives Ai[u>(i)] = X°D 1 / 3 [w(t)]/T 2 / 3 . Using the ex- 
pression for D{a) and averaging over w, we find 

fc.C, f» GS A>-^)(J^) 1/S >4 (6) 

For A7 1 [7i'(i)] ~ r the above procedure gives the order of 
magnitude estimate Ai ~ (X\(w)). Using the results on 
Kraichnan model, this gives Air ~ 1 at Gt ~ 1. Moder- 
ate order moments (w -1 / 3 ) and J °° S2(r)dr/r 2 entering 
G are expected to be well described by K41 at realistic 
Re. Indeed using for the estimates the multifractal model 
[2| to analyze J S2(r)dr /r 2 and express the anomalous 
exponents of w via empirical values of spatial anomalous 
exponents :18], one finds that intermittency becomes im- 
portant at Re well above 10 15 . Thus the above use of 
K41 for order of magnitude estimates is self-consistent. 
For G we obtain Gt - St 1/6 . We find the behavior of Ai 
similar to Kraichnan model. The region of weak inertia, 
Ai ~ Af lrt -C 1/t, is separated from the region of strong 
inertia Air ^> 1 where Eq. ((6]) holds, by the long crossover 
St > St 1 / 6 ~ 1 where Air ~ 1 is a slowly varying func- 
tion of St (note Air - ((D(u;)r) 1 /3A 1 [((Z)'( w ) T )-i/3)]) at 

\/St > 1). The limit St 1/6 > 1 at r < t L means very 
large Re and is of theoretical value mainly. Summarizing: 

Ai/Ai urt ~ St" 5/6 for St > 1, (7) 

The decay of Ai/A' ur6 at increasing r is faster than in 
Kraichnan model because of D — r -1 / 2 . 

We now consider the growth of R in the inertial range, 
at 77 < i? <C Zr- In contrast to Richardson's law for fluid 
particles R(t) ~ e 1/2 t 3 / 2 @,[Ti|, K41 dimensional analy- 
sis does not fix the separation law for inertial particles, 
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due to the additional time-scale r. We shall assume mod- 
erate i?(0) not to have mere ballistic motion, e. g. the 
analysis below applies to R(0) ~ r), R(0) ~ Ai?7, holding 
after separation at R <c r\. As we will see, i?(i) reaches 
i T within t ~ r, so to study separation at i?(f) <C Z T we 
assume £ <C r. The "friction" term fl/r in Eq. © pro- 
duces negligible effect over i -C r and can be omitted. 
Also £«(#) « + wt + R(t), t) - u(q(t) + wt, t), 

where w = w(Q). The correlation time t c (R) of Su(R) is 
due to the drift, t c (R) ~ < 7;^. As we verify later, 
the time-scale t c (R) of variations of R obeys t c (R) 3> 
t c (R). Proceeding like in the viscous range, we intro- 
duce Langevin description of Sui(R), substituting it by 
white noise Dij(R)jj, where (7i(i)7j(*0) = <>ij5(t — t'). 
Here D lk {R)D jk {R) = / dt([u t {R) ~ u t (Q)][u 3 {wt + 
R) — Uj(wt)}) to provide the correct dispersion of the 
time-averaged Su(R) [l7j]. We assumed for simplicity 
R/w < t R or (Ir/R) 1 / 3 > 1 (we use K41 as at R < if), 
so that the time correlations of 5u(R) are determined 
by the drift (cf. to r\ja < t v at i? < 77). The 
above Kraichnan model for particles is not the same as 
used usually to model turbulence in the inertial range 
t c (R) depends on R differently than t R . Noting 
D ik {R)D jk (R) ~ S 2 {R)t c {R) ~ S 2 (R)R/w, we conclude 
that the dependence on e and r in Ri = Dij(R,w)~/j/T 
is via single parameter e 2 / 3 /wt 2 ~ l¥ 3 /r 3 . Now di- 
mensional analysis is enough to fix the answer. We 
find t c (R) ~ r(i?// T ) 1 / 9 , so the applicability condition 
r c (i?) > iJ/iu gives (R/l T ) 8 / 9 < 1, close to just i? < / T . 
At rj ^> r c [i?(0)] the initial condition is forgotten (we 
assume explosive separation characteristic of the inertial 
range and R(t) depends only on t and l^ 3 /r 3 giving 

R(t) ~ l T (t/ T ) 9 , r c [i?(0)]<<<T. (8) 

In Kraichnan model, the power-law exponent for dis- 
persion of fluid particles grows indefinitely as the flow 
becomes less rough (for smooth flow separation is ex- 
ponential) [l9[. Thus in Eq. ((8|) exponent larger than 
in Richardson's law can be attributed to effectively 
smoother turbulence felt by particles. As R(0) > r], ob- 
servability of the power-law entails r 3> r c [i?(0)] > T c (rj). 
This gives (/ r /r 1 ) 1 /9 ^ St 1/6 > 1, equivalent to the nat- 
ural "forgetting" condition Af 1 < r. At St 1/6 ~ 1, the 
time of forgetting of the initial condition obeys t c [R(0)} ~ 
r so i?(t) at t -c t depends on the details of initial con- 
ditions. Eq. ([8|) then can be used as order of magnitude 
estimate at t ~ r giving -R(t) ~ / T . This is expectable - 
for fluid particles the time of separation to l T is of order r 
and determined by the stage of evolution with R{t) ~ l T 
where fluid and inertial particles behave similarly. 

To summarize, we have shown that the difference of 
the single particle velocity and the local velocity of the 
flow grows with inertia as the Lagrangian velocity incre- 
ment of turbulence at time r. While at St <§; 1 parti- 
cles disperse like fluid parfticles, at St 3> 1 there is a 



scale l T ~ r]St 3 ^ 2 below which particle dispersion obeys 
laws special for inertial particles and contribution of not 
too intermittent events can be described by white noise. 
The Lyapunov exponent Ai ~ A^" r6 /St 5 / 6 is estimated as 
the average of the Lyapunov exponent for constant drift 
problem over the drift velocity, the estimate becoming ex- 
act at St 1/6 > 1. In the inertial range, rj < R <C l T i the 
analogue of Richardson's law is R(t) ~ l T (t/r) 9 . The law, 
observable only at St 1 ^ 6 » 1, shows property expected 
at any St ^> 1: explosive separation to l T within t ~ r, 
closer than Richardson's law to exponential separation 
for smooth flows. The smoothing is by the effective time- 
averaging of turbulent velocity difference Su(R) driving 
the separation. Our treatment can be easily general- 
ized to incorporate constant gravitational or electric field, 
where constant drift problem emerges 0,113 (e.g. for 
gravity acceleration g one finds a — gr and Ai = Ai [gr]). 
The main qualitative result of the work is that collective 
drift of inertial particles through the flow makes their 
relative motion subject to Langevin description. 
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